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Abstract. The aim of this talk is to present the most recent advances 
in establishing plausible planetary system architectures determined by 
the gravitational tidal interactions between the planets and the disc 
in which they are embedded during the early epoch of planetary sys- 
tem formation. We concentrate on a very well defined and intensively 
i-C studied process of the disc-planet interaction leading to the planet mi- 

gration. We focus on the dynamics of the systems in which low-mass 
q planets are present. Particular attention is devoted to investigation 

t-H °f the role of resonant configurations. Our studies, apart from being 

ryj complementary to the fast progress occurring just now in observing the 

whole variety of planetary systems and uncovering their structure and 
origin, can also constitute a valuable contribution in support of the 
I missions planned to enhance the number of detected multiple systems. 
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jn 1 Diversity of planetary systems and migration-induced resonances 

1 I Multiple-planet systems are of particular interest to test theoretical models of 
planet formation and the evolution of planetary systems. Just a few days before 
this conference, Jason Wright (Wright et al. 2009) sent a paper in which he and 
his colleagues generated a catalogue of the 27 published multiple-planet systems 
around stars within 200 pc from the Sun to the other participants. The authors 
made several interesting observations about the systems' statistical properties. At 
least five of these systems harbour planets which are in or near mean-motion res- 
onance. This fact makes the importance of orbital migration even stronger than 
immediately after the discovery of giant planets very close to their host stars (e.g. 
Kley 2000; Papaloizou et al. I2007[l . These arguments, together with the increas- 
ing number of known planets with masses in the range of a few Earth masses, 
have stimulated our study of migration-induced resonances (e.g. Papaloizou & 
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Fig. 1. The initial configuration: two planets, with masses mi and mi respectively, in 
circular orbits around a central star with mass M, at distances r\ and r2, are embedded 
in a gaseous disc. 



Szuszkiewicz 2005, Podlewska & Szuszkiewicz I2008|) in systems containing low- 
mass planets. 

Our present goal is to investigate migration-induced architectures of planetary 
systems, focusing mainly on resonant configurations. In order to achieve this goal 
we have performed a series of numerical simulations of planets migrating in res- 
onance. The main tools used in our studies are two-dimensional hydrodynamic 
simulations, simple analytic modelling and N-body investigations. By construct- 
ing a simple analytic model, we were able to verify the reliability of our numerical 
calculations. By combining the hydrodynamic simulations with the N-body tech- 
nique we were able to follow the dynamical evolution of the planets for a substantial 
amount of time, comparable with the estimated lifetimes of gaseous discs. 



2 Initial configuration and computational set-up 

Our initial set-up shown in Fig. [T] includes the central star with a mass M + and two 
orbiting planets with masses vn\ and mi respectively. The two planets are embed- 
ded in a disc which is the source of planet orbit migration. They are initialised on 
circular orbits around a central mass which has a fixed value of 1 solar mass. The 
gravitational potential was softened with softening parameter b = 0.8H, where H 
is the semi-thickness of the disc. This results in the formation of an equilibrium 
atmosphere around the embedded planet which then does not accrete. This soften- 
ing also allows an adequate representation of type I migration in two dimensional 
discs (see e.g. Nelson & Papaloizou 12004)1 . The disc in which planets are ini- 
tially embedded has a surface density profile S(r) essentially flat (see Papaloizou 
&: Szuszkiewicz 120051 for a full description) . The planets are located in the flat 
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part of this distribution. We use four different values for the maximum value of 
the surface density So, namely the standard value of 2 x 10 3 kg/m 2 , attributed to 
the minimum mass solar nebula at 5.2AU which we denote Si, then So. 5 = 0.5Si, 
S 2 = 2Si and finally S 4 = 4S X . The radial boundaries were taken to be open. 

The calculations presented here were performed with a two-dimensional version 
of the Eulerian hydrodynamic code NIRVANA. For details of the numerical scheme 
and code adopted, see Nelson et al. (2000 ). 



3 A Super-Earth (mi = 4M®) and an Earth analogue (7712 = 1M©) 

Let us consider here a system of two planets with disparate masses, namely a 
Super-Earth with a mass of 4M® and an Earth analogue (1M®). The time-scale 
of inward migration for a low-mass planet on a circular orbit embedded in a disc 
with constant surface density can be approximated by the formula given by Tanaka 
et al. PTMII 



M, M, ■ 1 1 



^^brr MV 1 (3-D 



Here M. t is the mass of the central star, m p is the mass of the planet orbiting at 
distance r = r p , S p is the disc surface density at r = r p , c is the local sound speed 
and fip is the angular velocity at r = r p . The numerical coefficient W m is 0.3704. 

The density waves excited by a low mass planet with small eccentricity in the 
disc lead to orbital circularization (e.g. Artymowicz 1993; Papaloizou & Larwood 
2000) at a rate that can be estimated as (Tanaka & Ward 12004)1 



Tr 



W c \r p n p 



(3.2) 



Here, the numerical coefficient W c is equal to 0.289. 



It is expected from equation (3.1) that two planets with different masses will 
migrate at different rates. In our case, the Super-Earth in the external orbit 
will migrate faster than the Earth analogue. This has the consequence that their 
period ratio will evolve with time and may accordingly attain and become locked 
in a mean-motion resonance (Nelson & Papaloizou 12002) Kley et al. 2000) . We 
have considered the pair of planets to be evolving in discs with S = So. 5, S = 
Si, So — S2 and So — S4. The evolution was followed until a resonance was 
established. The results arc summarised in Fig. [2] (left panel) where the evolution 
of the ratio of the semi-major axes, which starts from the value 1.2 for all four 
cases, is shown. The fastest migration (steepest slope) corresponds to the case 
where the two planets are embedded in a disc with S = S4 and the slowest to the 
case of a disc with So = So. 5. The planets in the disc with S = S4 become trapped 
inap+l:p = 8:7 resonance. If the disc surface density is two times smaller, then 
a 7:6 resonance is attained. If it is four or eight times smaller, then the attained 
resonance is 5:4. These results are fully consistent with the idea that higher p 
resonances are associated with faster relative migration rates. Let us concentrate 
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Fig. 2. (left) The evolution of the ratio of semi-major axes for the two planets with 
masses mi = 4M® and ni2 = 1M®. Starting from the lower curve and going upwards, 
the curves correspond to the initial surface density scalings Eo = £4, So = £2, £0 = Si, 
and Eo = So. 5 respectively, (right) The semi-major axis ratio when mi = 7712 = 4M® 
with Eo = Ei (uppermost curve) and with Eo = E4 (two lower curves) 



on the 8:7 resonance, illustrating our combined methodology of hydrodynamic 
calculations (an explicit evolution of a 2D gaseous disc with embedded planets is 
calculated), simple analytic modelling and N-body techniques (the migration and 
eccentricity damping due to the gaseous disc is included in the form of analytic 
prescriptions) . 



3.1 Two-dimensional hydrodynamic simulations 

The whole evolution for the fastest migration rate is illustrated in Fig. [3] This 
shows the evolution of the individual semi-major axes, eccentricities, angle between 
apsidal lines and one of the resonant angles. The inner planet did not migrate 
significantly until an 8:7 resonance was attained at about 6000 years. Subsequently, 
the two planets migrated inwards together, maintaining the commensurability. 
During 18000 years of evolution the outer planet changed its location from a 
dimcnsionlcss radius of 1.2 to 0.9. The eccentricity of the outer planet increased 
slightly and that of the inner planet substantially, reaching an equilibrium value of 
0.04 at around 10000 years and at later times oscillating around this value. At the 
end of the simulation shown here, the ratio of eccentricities e\/&2 is equal to 0.125. 
The local peaks in the values of inner planet eccentricity occurring at 2000, 4000 
and around 5400 years correspond to the planet passing through the 5:4, 6:5 and 
7:6 resonances respectively. After about 8000 years the angle between the apsidal 
lines oscillates around 212° and the resonant angle around 264°. The amplitude 
of the oscillations for both angles decreases with time. 
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Fig. 3. The evolution of semi-major axes, eccentricities, angle between apsidal lines 
and resonant angle for two planets with masses, mi = 4M$ and m-2 = 1M® migrating 
towards a central star embedded in the disc with £ = £4, obtained by hydrodynamic 
simulations and N-body simulations (dotted lines in two upper panels and dark lines in 
two lower panels). 



3.2 Comparison with the analytic model 

It is of interest to compare the eccentricities obtained above when the planets are 
trapped in a commensurability with what is expected when resonant effects and 
disc tides are in balance. Papaloizou & Szuszkiewicz (2005) have given a simple 
approximate analytic solution for two migrating planets locked in a p + 1 : p 
commensurability In particular, they have shown that if eccentricities are not too 
large, then when they stop growing, they must satisfy 

i + ^_/i_e|\ /= M_ n/ (33) 

r cl t c2 min 2 a2 \t c l <^c2 J \ r migl t mig2 J o 

where / = TO2GSi/((p + l)(m2tli + 771^2))- Here the semi-major axes and eccen- 
tricities of the two planets (i = 1,2) are and e,. The migration rates (as- 
sumed directed inwards) and circularization times induced by the disc tides are 
tmigi = Wi / hi\ = |2aj/(3<x,)| = 2r r /3, and t c i = |e</e<| respectively. The mean 
motions are rij. 

In the simulations presented here e 2 » &\- We can simplify matters even 
further, setting e\ = 0. We can also assume that p is large, because in the present 
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example its value is 7. We also note that from equations (3.1l and (3.2) it turns 
out that t c — (T r /W c )(H/r) 2 . Thus we obtain 



el 



mi 
m 2 



1/2 



rn i 



0.578(p+1)(to 2 + toi) 



(H/rf 



(3.4) 



If we apply the above relation to the case illustrated in Fig. [3] for the two planets 
with masses m\ — 4M^ and mi = M® embedded in a disc with the aspect ratio 
H/r — 0.05 and initial surface density scaling Eo = £4, we obtain e 2 = 0.037 in 
reasonable agreement with the simulations. 



3.3 Comparison between hydrodynamic simulations and N-body calculations 

The hydrodynamic calculations discussed in the previous sections allow us to follow 
the migrations of planets for almost 2xl0 4 years. As a result, we have been able to 
simulate planets becoming trapped in resonances. The behaviour and stability of 
the resonance trapping vary with the planet masses and the surface density of the 
disc in which they are embedded. The outcomes of these simulations could be well 
matched to those of N-body integrations where we incorporate simple prescriptions 



for the migration and eccentricity damping given through equations (3.1 3.2 ). 
Using these expressions in the N-body code we have extended the hydrodynamic 
calculations for a longer period of time and studied the long term stability of 
the resonances that we found. As a first step we have adjusted the numerical 
coefficients in equations ( |3.1|3.2 1 in such a way that the hydrodynamic and N- 



body approaches give the same qualitative evolution. 

As an example, in Fig. [3] we show the results of the comparison for the case of 
two planets with masses IM9 and 4M^ respectively embedded in a disc with initial 
surface density scaling parameter So = S4. The numerical coefficients adopted 
were W m — 0.3647 and W c — 0.225. These results show good agreement with 
the hydrodynamic results and display the same 8:7 commensurability. The fitted 
coefficients were also reasonably close to those expected from analytic disc-planet 
interaction theory. 

The evolution was followed using the N-body approach for an additional 1.2 x 10 6 
years. The two planets remain in the 8:7 resonance during this time and there 
is no indication of any significant change in the monitored quantities for the last 
9 xlO 5 years. The long term evolution of this system is shown in Fig. [4] It is 
interesting to note that the equilibrium value of eccentricity of the inner planet is 
around e 2 = 0.03, in accordance with the value predicted by the simple analytic 
model discussed in Section 13.21 



4 Two Super-Earths (mi = m,2 = 4M e ) 

The resonant interaction must balance the tendency towards relative migration of 
the two planets. This is smaller when the planets have the same mass. In order 
to investigate this fact, two planets of mass 4M e were initiated close to a 3:2 
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time Tyearsl time [years! 

Fig. 4. The evolution of the semi-major axes and eccentricities for two planets with 
masses mi = 4Mg and mi = lMj migrating towards a central star embedded in the 
disc with E = E4 obtained from N-body simulations. 



resonance with r\ = 1.32 and r 2 = 1 in a disc with E = Si. The evolution of the 
semi-major axis ratio for the planet orbits is shown in Fig. [2] (right panel). We 
also performed simulations in a disc with Eo = £4 for the same masses starting 
at n = 1.23 and r 2 = 1.00 and at 7*1 = 1.20, r 2 = 1.00. The planets become 
trapped in the nearest available resonance which is a good indication of stability. 
The evolution of the equal mass pair of planets is discussed more extensively in 
Papaloizou & Szuszkiewicz (200.5)) and in Lacny & Szuszkiewicz (this volume). 



5 A gas giant and a Super-Earth (m\ = 333M®, m2 = 5.5M®) 

In the previous sections we have investigated systems in which both planets have a 
low mass. Here we extend our study to the case where one planet is a Super-Earth 
and the other a gas giant. The migration rates for different planet masses has been 
estimated by a number of authors, see the review by Papaloizou et al. (|2007|) . 
Their results are illustrated in Fig. [5] (left panel), where we plot the migration 
time of a planet as a function of its mass. There are two mass regimes which are 
of interest here, namely (0.1-30M®) and (150-1500MQ), for which in a typical 
protoplanetary disc we can talk about two different types of migration, called type 
I and type II respectively. The migration time for low mass planets embedded in 
a gaseous disc (type I migration) has been discussed above in Section [3] Type II 
migrators open a gap in the disc and their evolution is determined by the radial 
velocity drift in the disc v r . The migration time can be estimated as follows (Lin 
and Papaloizou 11993)1 

r 9r 2 

ru= r f = -^- (5.1) 

where v is a kinematic viscosity parameter. This has been shown in Fig. [5] (left 
panel) for r p = 5.2AU and different values of v (10~ 5 , 2 x 10~ 5 , 3 x 10"%..., 9 
x 10~ 5 and 10~ 6 ) expressed in dimensionless units. In drawing lines for a given 
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Fig. 5. (left) Comparison between the times of migration for planets with different masses 
located at 5.2 AU, embedded in a disc with a flat surface density distribution, the stan- 
dard value attributed to the minimum mass solar nebula at 5.2 AU, namely £ = 2000 
kg/m 2 and H/r = 0.05. (right) The evolution of the ratio of semi-major axes for the 
Jupiter and Super-Earth embedded in a gaseous disc. In the case of the upper curve the 
initial semi-major axis ratio is 1.45 and planets became locked into a 3:2 resonance. For 
the lower curve, the initial semi-major axis ratio is 1.24 and the attained commensura- 
bility is 4:3. The solid horizontal lines are the 3:2, 4:3 and 5:4 resonance widths in the 
case of circular orbits in the restricted three body problem (Wisdom, 119801 Lecar et al. 
12001"]) . 



viscosity parameter we have taken into account the condition for a gap opening in 
the disc which reads ^f- > ^f-- From this expression, it is clear that the bigger 

the viscosity v, the bigger the mass of the planet must be to be able to open a 
gap. We have chosen the mass of the type I and type II migrators to be 5.5M e 
and one Jupiter mass respectively. Their locations have been marked in Fig. [5] 
according to their masses and, in the case of the Jupiter, also according to the 
viscosity adopted in the disc (v = 10 -5 ). It is obvious from this figure that we 
should expect a convergent migration if the Jovian type planet is in the external 
orbit and the Super-Earth in the internal one. Starting from the configuration 
illustrated in Fig. [5] (left panel), we have studied the possible resonances in this 
system. In order to investigate how the outcome of such evolution depends on 
the initial planet configurations, we have performed our simulations with a wide 
range of planetary separations (see Table 1 in Podlewska & Szuszkiewicz 2008). 
The 2:1 commensurability was not taken into account because of computational 
time constraints. The resonances located closer to the Jupiter, such as 5:4, 6:5 or 
higher values of p, are not possible, because for such small separations the system 
becomes unstable and the Super-Earth is scattered from the disc. We have found 
that two outcomes are possible, namely either the Super-Earth is ejected from 
the disc or the planets become locked in 3:2 or 4:3 mean-motion resonances. In 
Fig. ^ (right panel) we show two examples of resonant trappings occurring when 
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the relative planet separation is 0.45 (upper curve) and 0.24 (lower curve). In the 
first case, the differential migration brought planets into a 3:2 commensurability 
and in the second case into a 4:3. The semi-major axis ratio of planets librates 
exactly within the regions of the resonance width marked in Fig. [5] (right panel) 
by solid, horizontal lines. Both commensurabilities shown in Fig. [5] are stable 
for the whole time of our simulations, namely 2500 orbits (the unit of time is 
defined as the orbital period on the initial orbit of the Super-Earth divided by 
2w). For full discussion of these results see Podlewska & Szuszkiewicz (J20J08J) and 
also Podlewska (this volume). 

6 Conclusions 

Particular attention has been devoted here to the investigation of the occurrence of 
resonant configurations in systems containing low-mass planets. Studies of com- 
mensurabilities have already been applied successfully to analyse the motion of 
pulsar planets and they proved to be a powerful tool in that context. They have 
the potential to be similarly useful in our search for Earth-like planets in other sys- 
tems. Our prediction of the occurrence of the commensurabilities in such systems 
may turn out to be particularly valuable in the case of the TTV (Transit Tim- 
ing Variation) observations, because the differences in the time intervals between 
successive transits, caused by planet-planet interactions, are largest near mean- 
motion resonances (e.g. Agol et al. 2005). Finally, detection of such resonances 
can also yield useful information about orbital migration as a process operating 
during planet formation. 
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